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Abstract 

We present constraints on the cosmological and reionization parameters based on 
the cumulative mass function of the Ly-o systems. We evaluate the formation rate 
of bound objects and their cumulative mass function for a class of flat cosmological 
models with cold dark matter plus cosmological constant or dark energy with con- 
stant equation of state encompassing different reionization scenarios and compare it 
with the cumulative mass function obtained from the Ly-a transmitted flux power 
spectrum. 

We find that the analysis of the cumulative mass function of the Ly-a systems 
indicates a reionization redshift Zr = 24.2 it 4 (68%CL) in agreement with the value 
found on the basis of the WMAP anisotropy measurements, setting constraints on 
the amplitude of the density contrast, as = 0.91 ± 0.04 (68%CL), similar to those 
derived from the X-ray cluster temperature function. 

Our joint analysis of Ly-a cumulative mass function and WMAP anisotropy mea- 
surements shows that the possible current identification of a running of the slope, 
dns/dlnk ^ 0, at A;p=0.05Mpc'~-'^ (multipole I 700) is mainly an effect of the exist- 
ing degeneracy in the amplitude-slope plane at this scale, the result being consistent 
with the absence of running, the other constraints based on WMAP data remaining 
substantially unchanged. 

Finally, we evaluate the progress on the determination the considered parameters 
achievable by using the final temperature anisotropy data from WMAP and from 
the forthcoming Planck satellite that will significantly improve the sensitivity and 
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reliability of these results. 

This work has been done in the framework of the Planck LFI activities. 

Key words: Cosmology: cosmic microwave background - large scale structure - 

dark mattor 
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1 Introduction 



The cosmological density background is assumed to have been seeded at some 
early epoch in the evolution of the Universe, inflation being the most popular 
of the current theories for the origin of the cosmological structures (see e.g. 
Kolb & Turner 1990, Linde 1990). 

Measurements of the Cosmic Microwave Background (CMB) anisotropies present 
at the epoch of the recombination at large scales {R^8h~^Mpc), always 
outside the horizon during the radiation-dominated era, carry precious in- 
formation on the power spectrum of density field. Moving to intermediate and 
small scales the background density field encodes information related to the 
non-linear evolution of the gravitational clustering, the time evolution of the 
galaxy bias relative to the underlying mass distribution (see e.g. Hoekstra et 
al. 2002, Verde et al. 2002), and the magnitude of the pecuhar motions and 
bulk flows in the redshift space (Kaiser 1987). At these scales the spectrum 
of the background density field can be constrained by complementing CMB 
measurements with other astronomical data (see Kashhnsky (1998) for the 
evaluation of the spectrum of the density field from a variety of astronomical 
data on scales l/i^^Mpc < R < lOO/i^^Mpc ). 

In particular, on scales less then 5 h~^Mpc the spectrum of the density field 
can be constrained by using observations of the spatial distribution of high- 
redshift collapsed objects as quasars and galaxies and information on their 
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macroscopic properties (Efstathiou & Ress 1988, Kashlinsky & Jones 1991, 
Kashlinsky 1993). The detection of high-z quasars by the Sloan Digital Sky 
Survey (Becker et al. 2001, White et al. 2003, Fan et al. 2003) and by the Keck 
telescope (Vogt et al. 1994, Songaila & Cowie 2002) as well as the detection 
of high-redshift galaxies (see Kashhnsky (1998) and the references therein) 
are indications about the existence of such early collapsed objects at redshifts 
between 2 and 6, on mass scales of ~ lO^^M©. 

The rms mass fluctuations over a sphere of radius of 8/i~^Mpc, erg, fixes the 
amplitude of the density field power spectrum and then determines the red- 
shift of collapse, Zc, of an object with a given mass-scale. Consequently, the 
observation of such objects can set significant constraints on different cosmo- 
logical theories. 

Recently, the WMAP ^ team (Spergel et al. 2003, Verde et al. 2003) high- 
lighted the relevance of complementing the CMB and LSS measurements with 
the Ly-a forest data (the absorptions observed in quasar spectra by the neu- 
tral hydrogen in the intergalactic medium) in constraining the cosmological 
parameters (Croft et al. 1998; Zaldarriaga et al. 2003 and references therein) 
and the shape and amplitude of the primordial density field at small scales. 

In this paper we investigate the dynamical effects of the non-linear evolution 
of the density field imphed by the high-z collapsed objects as diagnostic of the 
reionization history of the Universe. The reionization is assumed to be caused 

^ http: / /lambda.gsfc. nasa.gov 
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by the ionizing photons produced in star-forming galaxies and quasars when 
the cosmological gas falls into the potential wells caused by the cold dark mat- 
ter halos. In this picture the reionization history of the Universe is a complex 
process that depends on the evolution of the background density field and of 
the gas properties in the intergalactic medium (IGM) and on their feedback 
relation. The evolution of the background density field, that determines the 
formation rate of the bound objects, is a function of the grow rate of the 
density perturbations and depends on the assumed underlying cosmological 
model. The evolution of the gas in the IGM is a complex function of the gas 
density distribution and of the gas density-temperature relation. The latter 
is related to the spectrum (amplitude and shape) of the ionizing radiation, 
to the reionization history parametrized by some reionization parameters (the 
reionization redshift, Zr, and the reionization temperature, T^), and to the as- 
sumed cosmological parameters. 

As shown by hydrodynamical simulations (Cen et al. 1994; Zhang et al. 1995; 
Hernquist et al. 1996; Theuns et al. 1998), the gas in the IGM is highly inho- 
mogeneous, leading to the non-linear collapse of the structures. In this process 
the gas is heated to its virial temperature. The photoionization heating and 
the expansion cooling cause the gas density and temperature to be tightly 
related. Finally, the temperature-mass relation for the gas in the IGM at the 
time of virilization determines the connection between the gas density and the 
matter density at the corresponding scales. 
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Taking advantage of the results of a number of hydrodynamical simulations 
(Cen et al. 1994; Miralda-Escude et al. 2000; Chiu, Fan & Ostriker 2003) 
and semi-analytical models (Gnedin & Hui 1998; Miralda-Escude, Haehnelt 
& Rees 2000, Chiu & Ostriker 2000), we re-assess in this paper the possibility 
to use the mass function of Ly-a systems to place constraints on spatially 
flat cosmological models with cold dark matter plus cosmological constant or 
dark energy, encompassing different reionization scenarios. We use the Press- 
Schechter theory (Press & Schechter 1974) to compute the comoving number 
density of Ly-a systems per unit redshift interval from the Ly-a transmit- 
ted flux power spectrum (Croft et al. 2002) and examine the constraints on 
the cosmological parameters. In our analysis we take into account the con- 
nection between the non-linear dynamics of the gravitational collapse and the 
properties of the gas in the IGM through the virial temperature-mass rela- 
tion. We address the question of the consistency of the WMAP and Ly-a 
constraints on the cosmological parameters, paying a particular attention to 
the degeneracy between the running of the effective spectral index and the 
power spectrum amplitude. Finally, we discuss the improvements achievable 
with the final WMAP temperature anisotropy data and, in particular, the 
impact of the next CMB temperature anisotropy measurements on-board the 
ESA Planck ^ satellite. 



^ http: / / astro.estec.esa.nl/Planck/ 
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2 Ecirly object formation mass function 



2.1 Formation rates 



The most accurate way used to assess the formation rates of the high-z col- 
lapsed objects is based on numerical simulations. A valid alternative is offered 
by the Press-Schechter theory (Press & Schechter 1974, Bond et al. 1991) ex- 
tensively tested by numerical simulations for both open and flat cosmologies 
(Lacey & Cole 1994, Eke, Cole & Frenk 1996, Viana & Liddle 1996). 

According to the Press-Schechter theory, the fraction of the mass residing in 
gravitationally bounded objects is given by: 



/coii(^) ~ \ -^^ — 7 / exp 

Oc(Z) 



d5,{z). (1) 



Here Sc{z) is the redshift dependent density threshold required for the collapse; 
Rf is the filtering scale associated with the mass scale M = iirRj^pb/S, pb being 
the comoving background density; ^{Rf, z) is the rms mass fluctuation within 
the radius Rf-. 

a'{Rf,z)^ l—A'{Kz)W'{kRf), (2) 



where W(x) is the window function chosen to fllter the density fleld. For a 
top-hat filtering W{x) = 3(sina; — xcosx)/x^ while for a Gaussian filtering 
W{x) = exp(— x^). As we do not find significant differences between the re- 



suits derived by adopting the two considered window functions for some rep- 
resentative cases, we present in this work the results obtained by using the 
Gaussian smoothing. In the above equation A^(/c,2;) is referred as the power 
variance and is related to the matter power spectrum P{k, z) through: 



Motivated in the framework of the spherical collapse model and calibrated by 
N-body numerical simulations, the linear density threshold of the collapse 5c{z) 
was found to vary at most by ~ 5 % with the background cosmology (see e.g. 
Lilje 1992, Lacey & Cole 1993, Eke, Cole & Frenk 1996). However, the choice 
of 5c{z) depends on the type of collapse. For the spherical collapse the standard 
choice is (5c(0) — 1.7±0.1 while for pancake formation or filament formation its 
value is significantly smaller (Monaco 1995). For the purpose of this work we 
assume that the collapse have occurred spherically and use 5c(0) = 1.686 that 
is the conventional choice for the case of the flat cosmological models (Eke, 
Cole & Frenk 1996). The time evolution of 5c depends on the background 
cosmology: 



where D[z) is the linear growth function of the density perturbation, given in 
general form by (Heath 1977, Carroll, Press & Turner 1992, Hamilton 2001): 



^\k,z)^^k'P{Kz). 



(3) 



D{z) 



(4) 



D(0)' 




(5) 
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Here a — {1 + z)~^ is the cosmological scale factor normalized to unity at 
the present time {gq = 1), flrn is the matter density energy parameter at the 
present time and /(a) specifies the time evolution of the scale factor for a 
given cosmological model: 



da Ho 
dt f{a) 



1 + 



l-de 



-1/2 



(6) 



In the above equation Hq is the present value of the Hubble paramenter, 
Q^de is the present value of the energy density parameter of the dark energy, 
w = p/ Pde ^ —1 defines the dark energy equation of state, and Vtm is the 
matter energy density parameter. Equation (6) reduces to that of a ACDM 
model for w — —1. 

The comoving number density of the gravitationally collapsed objects within 
the mass interval dM about M at a redshift z is given by (Viana & Liddle 
1996): 



n{M,z)dM = -1 



da{Rf, z) 



7rMa'^(Rf,z) dM 



exp 



2a^(Rf,z) 



dM. (7) 



We are interested in the formation rate of the high-z collapsed objects at 
a given redshift. According to Sasaki method (Sasaki 1994), the comoving 
number density of bounded objects with the mass in the range dM about 
M, which virilized in the redshift interval dz about z and survived until the 
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redshift Zf without merging with other systems is given by: 



N(M, z)dMdz 



n{M,z) da{Rf,z) 



a'^{Rf, z) cr(Rf, . 



dz 



a{Rf, Zf) 



dMdz, (8) 



where: 

The total comoving number density of the bounded objects per unit redshift 
interval with the mass exceeding M (the cumulative mass function) is given 
by: 



N(> M)^ J N(M, z)dM. 

M 



(9) 



2.2 Mass-temperature relation for the virilized gas in the IGM 

The fraction of the mass of the gas in collapsed virilized halos can be calculated 
if the probability distribution function (PDF) for the gas overdensity is known: 



oo 

/eoll(^)= / APv(A)d 



A. 



(10) 



Here -Py(A) is the volume-weighted PDF for the gas overdensity A = Pg/pbar, 
where pg is the gas density and pbar is the mean density of baryons. In the 
above equation Ac{z) is the halo density contrast at virilization. Based on hy- 
drodynamical simulations, Miralda-Escude et al. (2000) found for the volume- 
weighted probability distribution, Pv{A), the following fitting formula: 



Pv{A)dA = Aexp 



- Co) 
2(25o/3)2 



A-f^dA. 



11 



In this equation 5o is the hnear rms gas density fluctuation and /3 is a parameter 
that describe the gas density profile ( Pg ~ for an isotermal gas). As shown 
by the numerical simulations (see Table 1 from Chiu, Fan & Ostriker 2003) 
the redshift evolution of Sq depends on the underlying cosmological model. 
Assuming the same fraction of baryons and dark matter in collapsed objects, 
we compute the redshift dependence of So on the cosmological parameters by 
using an iterative procedure (Chiu, Fan & Ostriker 2003) requiring the equality 
of the equations (1) and (10) representing the collapsed mass fraction. The 
parameters A and Co were obtained by requiring the normalization to unity of 
the total volume and mass. For the redshift dependence of /? in the considered 
redshift range we take the values P min[2.5, 3.2 — 4.73/ (1 + z)] obtained by 
Chiu, Fan & Ostriker (2003) through a fit to their hydrodinamical simulations. 

As we are interested to apply equation (9) to the virilized gas, we need to know 

the temperature-density and the mass-temperature relations for the virilized 
gas in the IGM. 

The temperature-density relation is determined by the reionization scenario 
and the underlying cosmological model. Hydrodynamical simulations can pre- 
dict this relation accurately, but the limited computer resources restrict the 
number of cosmological models and reionization histories that can be studied. 
For this reason we evaluate the temperature-density relation at the redshifts 
of interest by using the semi-analytical model developed by Hui & Gnedin 
(1997) that permits to study the reionization models by varying the ampli- 
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tude, spectrum, the epoch of reionization and the underlying cosmological 
model. According to this model, for the case of uniform reionization models, 
the mean temperature-density relation is well approximated by a power-law 
equation of state that can be written as: 

T^Toil + Ay-\ (12) 



where To and 7 are analytically computed as functions of the reionization 
temperature T^, the reionization redshift Zr, the matter and baryon energy 
density parameters ftm and ^bar, and the Hubble parameter Hq. 

According to the virial theorem (Lahav et al. 1991, Lilje et al. 1992), the virial 
mass-temperature relation at any redshift z can be written as (Eke, Cole & 
Frenk 1996; Viana & Liddle 1996; Kitayama & Suto 1997; Wang & Steinhardt 
1988): 

-3/2 

,(13) 

where /c^ is the Boltzmann constant, T^j^ is the temperature of the virilized 
gas, Ac is the density contrast at virilization and fp — fulJ'/P, where /„ is the 
fudge factor (of order of unity) that allows for deviations from the simplis- 
tic spherical model and jj, is the proton molecular weight. Different analyses 
adopting similar mass-temperature relations disagree on the value of fp be- 
cause of the uncertainties in the numerical simulations. We adopt here = 1, 
as indicated by the most extensive simulation results obtained by Eke, Cole 
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0/ ^ 



W-^H-^Mq l 0.944keVy 



[(1 + ^)3noAe 



-1/2 



2^de(-g) 



& Prenk (1996). 

Throughout this paper we consider that the virihzation takes place at the 
collapse time, t{zc), that is half of the turn-around time: t{zc) — t{zta)/'^, 
Zta being the redshift at which R{zta) = (i? is the radius of a spherical 
overdensity) . 

For ACDM models the density contrast at virilization is a function of flrn 
only. For quintessence models the density contrast at virilization becomes a 
function of fl^ and w and can be written as (Wang & Steinhardt 1988): 



pb{Zc) \RvirJ \l + Zc 

where: Ci^zta) = Pciust{zta) / Ph{zta)] Rta and R^ir are the radius at Zta and Zc 
respectively: 

azta) = (37r/4)20^(z,„)-°-^9+°-26^'"(-)-°-°6-, (15) 

Ryir 1 Tjy 1 2 



(16) 



where 77t = 2C ^^de{zta)l^m{zta) and 77^ = 2C ^[(l+2;c)/(l+2;ta)]^Ode(2;c)/f^m(2;c)- 
The energy density parameters and the Hubble parameter evolve with the scale 
factor according to: 

n^{a)^^^, f),^(a) = ___, i/(a) = __, (17) 
where /(a) is given by the equation (6). 

The scale factors for collapse, Oc, and turn-around, ata, was computed from 
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the spherical collapse model (Lahav et al. 1991, Eke, Cole & Prenk 1996). 
For any region inside the radius R that was overdense by Aj with respect to 
the background at some initial time ti corresponding to the redshift Zi, we 
solve the set of equations: 



ata sta flc Sta 

j f{a)da= j g{s)ds, j f{a)da = j g{s)ds, (18) 







where: 



ds Hi 



and g{s) 



1 + Qiil + AO ( - - 1 ) + ^de,i{^ - 1) 



-1/2 



.(19) 



dt g{s) 

In the above equations, solved by using an iterative procedure, s — R/Ri 
is the scale factor of a spherical perturbation with the initial radius Ri and 
Sta — Rta/Ri is its scalc factor at the turn-around (see Appendix A in Eke, 
Cole & Frenk 1996). The average density perturbation inside the radius R is 
given by: 

A(i?, z) = ^ j R:'5{R, z)d R, (20) 



where d{R, z) is related to the power spectrum of the density field P{k, z) — 
|5fc(2;)P through: 

^) = / 5,(^)e-'=«d^ k. (21) 



As initial conditions for the spherical infall we choose the epoch given by 
Zi — 1100 when the growth of perturbations is fully determined by the linear 
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theory. The initial values of the parameters Hi — H{zi), Qj = ^{zi) and 
flde,i = ^de{zi) are given by the equation (17) and for the normalization of the 
density field at Zi, as{zi), we take: 

We adopt for erg at the present time the value obtained from the analysis of the 
local X-ray temperature function for flat cosmological models with a mixture 
of cold dark matter and cosmological constant or dark energy with constant 
equation of state (Wang & Steinhardt 1998): 

as = (0.50 - 0.ie)Q-'^(^'"'®) , (22) 

where: 

-f{nm, e) = 0.21 - o.22w + o.ssnm + 0.259 

and 

e = (n5-l) + (/i-0.65). 
3 Results 

3. 1 Cosmological constraints from Ly-a observations 



The study Ly-a transmitted flux power spectrum has become increasingly im- 
portant for cosmology as it is probing the absorptions produced by the low 
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density gas in voids or mildly overdense regions. This gas represents an ac- 
curate tracer of the distribution of the dark matter at the early stages of the 
structure formation. One of the most important application is to recover the 
linear matter power spectrum Phik) from the flux power spectrum Ppik) and 
inferring the cosmological parameters of the underlying cosmological model. 
On the observational side there are recent analyses by McDonald et al. (2000) 
and Croft et al. (2002) that obtain results for the transmitted flux power 
spectrum Ppik) in agreement with each other within the error bars. Two dif- 
ferent methods have been proposed to constrain the cosmological parameters: 
McDonald et al. (2000) and Zaldarriaga et al. (2001) directly compare Ppik) 
with the predictions of the cosmological models, while Croft et al. (2002) and 
Gnedin & Hamilton (2002) use an analytical fitting function to recover the 
matter power spectrum, Pl(^), from the flux power spectrum Ppik). The 
WMAP team (Verde et al. 2003) used the analytical fitting function obtained 
by Gnedin & Hamilton (2002) to convert Ppik) into PL{k). 
In a recent work, Seljak, McDonald & Makarov (2003) investigate the cos- 
mological implications of the conversion between the measured flux power 
spectrum and the matter power spectrum, pointing out several issues that 
lead to the expansion of the errors on the inferred cosmological parameters. 

We compute the total comoving number density, A^L^z-a, of Ly-o; systems per 
unit redshift interval that survived until z = 2.72 without merging with other 
systems from the flux transmission power spectrum, Ppik), obtained by Croft 
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at al. (2002) for their fiducial sample with mean absorption redshift z — 2.72. 
We then compare NLy_a with the theoretical predictions for the same func- 
tion, Nth, obtained for a class of cosmological models encompassing the dark 
energy contribution with constant equation of state and different reionization 
scenarios. 

Our fiducial background cosmology is described by a fiat cosmological model, 
VLq = il^ + fide = 1) with the following parameters at the present time: 
Qbarh'^ = 0.024 ± 0.001, Qrrih^ = 0.14 ± 0.02, h = 0.72 ± 0.05, as indicated 
by the best fit of the power law ACDM model of WMAP data (Spergel et al. 
2003). In our analysis we allow to vary the primordial scalar spectral index 
Us, the parameter w describing the equation of state for the dark energy, the 
reionization redshift Zr, and the reionization temperature T^. We assume adi- 
abatic initial conditions and neglect the contribution of the tensorial modes. 



Our parameter vector p = (n^, w, Zr, T^.) has four dimensions. We create a grid 
of model predictions for the each choice of the parameters in the grid: 

• (0.7, 0.8, 0.85, 0.99, 1.1, 1.15, 1.2, 1.3) 

• (-0.4, -0.5, -0.6, -0.65, -0.7, -0.75, -0.8, -0.85, -0.9, -0.95, -1) 

• Zr^ (5, 6, 8, 10, 12, 14, 16, 18, 20, 22, 24, 26, 28, 30, 35, 40) 

• Tr,A = (1.6, 1.8, 1.9, 2, 2.1, 2.2, 2.3, 2.4, 2.5, 2.6, 2.8, 3) 

Here Tr^^ is the reionization temperature in units of lO^K. 
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Fig. 1. The dependence of the virihzed gas properties on the cosmological and 
reionization parameters: Ug = 0.99, Zr = 10, w = —1, T^^i = 2.5 (solid Unes), 
Ug = 0.99, Zr = n, w = —1, Tr^4 = 2.5 (dashed hnes), Ug = 0.99, Zy = 10, w = — 1, 
Tr,i = 2 (dot-dashed Unes), Ug = 1.1, Zr = 10, w = —1, Tr^4 = 2.5 (small dot-dashed 
lines), Hg = 1.1, Zr = 10, w = —0.7, Tr^4 = 2.5 (dotted lines). rr,4 represents the 
reionization temperature in units of lO^K. See also the text. 
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The density perturbations Sk{zi) at the initial redshift Zi — 1100 and the mat- 
ter transfer function at z = 2.72 was computed for each set of parameters 
in the grid by using the CMBFAST code version 4.2 (Seljak & Zaldarriaga 
1996). Then we evaluate the linear matter power spectrum, Phik), with the 
appropriate normalization. 

The Ly-a transmission power spectrum for the fiducial sample at ^ = 2.72 
proves linear scales in the range R = (11.4 — 1750) km/s (see Table 2 from 
Croft et al. 2002). For the purpose of this work we consider scales up to 
791 km/s, which are non-linear today. We compute the averaged density per- 
turbation Aj inside each scale R at the initial redshift Zi and the correspond- 
ing redshift of collapse, Zc, as described in the previous section. Then, as- 
suming that the virilization takes place at the collapse time, we evaluate the 
temperature-density relation at Zc and compute the appropriate virial mass at 
z — 2.72 by using the mass-temperature relation. The virial mass obtained in 
this way is related to the filtering scale through M^ir = {4:7i/3)Rjpb{z) with 
Pb{z) — {3 Hq/ 87:0)0,0(1 + z)^. We found that the filtering scale obtained in 
this way, Rf pa Rj/2, is proportional to the Jeans scale, Rj, corresponding 
to the the Jeans mass Mj = 1.5T4(1 + z^y^/^n-^/HO^^h-^MQ (T4 being the 
gas temperature in units of 10^ K), but the exact relation depends on the 
cosmological model and reionization parameters. 

Figure 1 presents few dependences of the gas properties on the background 
cosmology and reionization parameters. Panels a) and b) show the dependence 
of the gas temperature Tq and of the parameter 7 on the redshift of collapse 
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Zc- Panel c) and d) present the dependence of the virihzed mass on the hnear 
scale R and on the mass scale M{R), as an indication on the fraction of the 
virilized mass at the given scale. In panel e) we show the dependence of the 
clumping factor C —< Pg > / < pbar >^ on the virial mass. Panel f) presents 
the dependence of the filtering scale Rf on the hnear scale R. 
We note that the filtering scale obtained in this way depends on our param- 
eter vector: Rf = Rf{k, p), where k is the wave number corresponding to the 
hnear scale R. 

For each choice of the parameters in our simulation grid we compute NLy^a ac- 
cording to the equations (7) - (9) by filtering the transmission power spectrum 
Ppik) at each wave number k with the corresponding filtering scale Rf{k, p) 
and apply the same procedure to compute Nth from the linear power spectrum 
PL{k). We compare N^y^^ and Nth by computing a Gaussian approximation 
of the likelihood function: 

17 

L{NLy-a;p) oc ]J = exp 
1=1 

where i runs over different points [we are using first 17 bins from the Ppik)] 
and (7i is the error bar associated to NLy^o. on each point. The full chi-squared 
goodness of fit is ~ — 21nL with 12 degrees of freedom (ndf). We define the 
confidence level (CL) as the upper tail probability of the chi-squared distri- 
bution and calculate x^(CL, ndf) for a given CL by inverting the chi-squared 
distribution. 
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(23) 



Fig. 2. The cumulative mass function of the Ly-a systems per unit redshift interval 
at z = 2.72 (filled and open circles) compared with the corresponding theoretical 
predictions (solid and dashed lines). See also the text. 
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Figure 2 presents N^y^a and Nth functions obtained for few choices of the 
parameters in the simulation grid. In each of the panels a), b), and c) we 
indicate in the top the values of the parameters in common to the two curves 
reported in each of the panels a), b), and c) are indicated in the top while 
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the values of the parameters associated only to one of the two curves of each 
of the panels a), b), and c) are indicated in the bottom together with the 
corresponding confidence level obtained from our analysis. Panel d) shows an 
example of two quasi-degenerated models and their predictions for NLy^o.- 

Each simulated -Pl(A;) z — 2.72 was parametrized by its effective slope 
(Peacock & Dodds 1996), the running of the slope and the power variance at 
the pivot point /c*=0.03 s/km: 

n*Ak;) = ^^^{k = k;/2), (24) 



{k = k;/2) , (25) 



Al{k;)^2n'k'P4k)U;. (26) 

Assuming full ionization (ionization fraction Xe—1) for an easier comparison 
with the results of the WMAP team, we compute for each model in our grid 
the optical depth to the last scattering, r, and the normalization of the matter 
power spectrum at the present time in terms of as- 

Once we have computed the value of the ^oi every choice of the parameters 
of our simulation grid we marginalize along one direction at a time to get the 
four- dimensional constraints on the considered parameters. 
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In Table 1 we report the constraints ^ we have obtained at 68% CL and 95% 
CL. Figure 3 presents our 68% CL contour in n*g - plane, indicating the 
best fit values and in Figure 4 we compare the best fit linear power spectrum, 
Pl(A;), at ^ = 2.72 (68% CL) with the transmission power spectrum Pf(/c). 
The error bars include the statistical error and a systematic error due to the 
undeterminations in ft^n, ^de and Hq. 

Our values for n*fj and are in a good agreement with those obtained by 
Seljak, McDonald & Makarov (2003) indicating the same degeneracy direction 
in n*g - A^ plane, but are only marginally consistent with the similar values 
obtained by Croft et al. (2002). 

We find that the analysis of the cumulative mass function of the Ly-a systems 
indicates a reionization redshift in agreement, within the error bars, with the 
value found on the basis of the WMAP anisotropy measurements, setting 
constraints on the amplitude of the density contrast similar to those derived 
from the X-ray cluster temperature function. 

3.2 Combined CMB and Ly-a analysis 

By jointly considering the WMAP anisotropy data and the Ly-a observations, 
we investigate the cosmological constraints obtained on the basis of the anal- 
^ We also verified that the final results do not change significantly by using the 
transmission power spectra from McDonald &: Miralda-Escude (1999) at z = 2.4,3 
and 3.9. 
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Table 1 

Cosmological parameter constraints from Ly-a cumulative mass function. The val- 
ues of n*g, dn*Q/d\nk and are evaluated at the pivot point A:*=0.03 s/km and 
z = 2.72; Tr^4 is the reionization temperature in units of lO^K. 



Ly-a 


Parameter 


68% CL 


95% CL 


ris 


1.001±0.034 


0.975±0.047 




2.317±0.205 


2.244±0.227 




24.195±3.976 


22.313±4.814 


W 


-0.689±0.141 


-0.729±0.144 


/To 


n qi 1 _|_n qqc 

\J .iJ J- J \—\J •\J*JLJ 


n qiq+n nQO 




0.148±0.035 


0.132±0.041 




0.786±0.041 


0.805±0.047 




-2.550±0.034 


-2.576±0.047 


(in*g / dink 


-0.017±0.004 


-0.018±0.005 




0.666±0.113 


0.618±0.131 



"Assumes ionization fraction, Xe=l. 
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Fig. 3. Constraints at 68% CL on the effective slope n*g and power variance 
at A;*=0.03 s/km and z = 2.11 from: Ly-a cumulative mass function analysis (thin 
solid contour and filled square), WMAP anisotropy measurements (dashed contour 
and filled circle) and the joint WMAP and Ly-a analysis (thick solid - green - 
contour and open circle). See also the discussion in Sect. 3.2. The similar values 
obtained by Croft et al. 2002 (filled triangle) are also indicated. 

1.2 
1 

0.8 

^ 0.6 
0.4 
0.2 

-2.7 -2.65 -2.6 -2.55 -2.5 -2.45 -2.4 -2.35 

ysis of the Ly-a cumulative mass function. To this purpose, we use here only 
the accurate WMAP temperature anisotropy (TT) power spectrum. In fact, 
although crucial to probe the cosmological reionization, the inclusion of the 
polarization (ET) power spectrum derived from WMAP [see also the DASI 
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Fig. 4. The best fit linear power spectrum Phik) at z = 2.72 obtained at 68% CL 
from the analysis of the Ly-a cumulative mass function (open circles and solid line) 
compared with the transmission power spectrum Ppik) from Croft et al. 2002 (filled 
circle) . 
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detection/upper limit on (E and B) polarization power spectrum; Kovac et 
al. 2002] does not change significantly our quantitative results, as we have 
verified for a representative set of cases. We will take into account the polar- 
ization information in future works. We ran the CMBFAST code v4.2 with 
the CO BE normalization option to generate the CMB temperature anisotropy 
power spectra for the considered grid of parameters and then renormalized 
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each computed power spectra, C^, of the grid to minimize the when com- 
pared to the WMAP data (we find that this renormahzation does not change 
appreciably the final best fit and error bar results). We vary n^, w and in 
the same range as in the previous analysis, but for this case we also allow to 
vary the effective running of the slope dnsld\vik\k^ at = 0.05Mpc~^, the 
same pivot wavenumber used in the analysis by the WMAP team (Spergel et 
al. 2003). 

Our parameter vector has also four dimensions: p = {ns,w, Zr,dns/dlnk), 
where dug/dhik was free to vary in the range [—0.1,0.1] with a step of 0.01. 
We compute the for each choice of the parameters in the grid, comparing 
the simulated CMB temperature anisotropy power spectra with the WMAP 
anisotropy power spectrum, following the same procedure as in the previous 
analysis. In addition, for each choice of parameters in the grid we generate 
the matter transfer function and evaluate the linear matter power spectrum 
Phik) at the present time with the normalization given by the equation (22). 
For each Phik) we compute the effective slope nes{kp) and the running of the 
slope dries/dlnk at the pivot wavenumber kp = 0.05Mpc~^, as given by equa- 
tions (24)-(26) by only replacing k* with kp. Note the difference between the 
two definitions of the running of the slope: dus/ dink is the matter power law 
free parameter (see Spergel et al. 2003 and the CMBFAST code v4.2) while 
dries I dh\k is obtained from the matter transfer function shape. 
We generate Monte Carlo Markov Chains (MCMC) using the combined 
obtained from the Ly-a cumulative mass function and WMAP anisotropy 
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Table 2 

Cosmological constraints from the joint WMAP and Ly-a analysis. The val- 
ues of Weff, dus/dlnk and dries/dlnk are evaluated at the pivot wavenumber 
A;p=0.05Mpc""^; Tr^4 is the reionization temperature in units of lO^K. 

WMAP+Ly-a 



Parameter 68% CL 95% CL 



ris 


1.034±0.043 


1.011±0.059 


Tr,4 


2.320±0.194 


2.228±0.229 




26.131±3.105 


25.495±3.728 


in 

UJ 




_r\ oqi _|_n 1 1 r 


0-8 


0.945±0.035 


0.949±0.032 




0.167±0.028 


0.159±0.032 




0.800±0.037 


0.810±0.040 


"eff 


-2.151±0.045 


-2.154±0.047 


dries 1 d Infe 


0.018±0.004 


0.017±0.005 


dnsjdhxk 


0.023±0.022 


0.024±0.023 



"Assumes ionization fraction, a;e=l. 
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Fig. 5. Panel a): constraints in UeS - dus/dlnk plane from MCMC with dng/cUnk ^ 0. 
Panels b) and c): constraints in neff - crge"^ plane and crge"^ - c/nefr/c?lnfc plane 
from MCMC with drieff /dlnA; = 0. All contours and error bars are at 68% CL; neff, 
dus/dlnk and dnes/dlnk are obtained at kp = 0.05Mpc~^. 

-1 

xlO 




-2.4 -2.2 -2 -2.4 -2.2 -2 0.6 0.8_ 1 

Heff ri,(, CTge ^ 

analysis (the number of model elements is of about 3 x 10^). In this way we 
sample the chi-squared probability distribution in the combined parameter 
space. We run the MCMC with dus/dlnk = and varying dUs/dlnk in the 
range mentioned above. 

Before of discussing the results obtained by using the pivot wavenumber 
A;p=0.05Mpc~^ we briefly report on the results we derived by using the pivot 
point A;;=0.03 s/km as in Sect. 3.1 but by exploiting only the WMAP (TT) 
power spectrum or WMAP combined to the Lj-a information. They are shown 
again in Figure 3: note how for this pivot point choice the poor sensitiv- 
ity of WMAP at the small scales accessible to Ly-a observations does not 
improve but slightly worses the parameter recovery, the two kinds of obser- 
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vations having similar degeneracy directions in the n*ff — plane. On the 
contrary, the situation improves by considering a pivot wavenumber (namely 
at A;p=0.05Mpc~^) at larger scales. We report in Table 2 the results obtained 
from the joint WMAP and Ly-a analysis from the MCMC with dng/ dhik ^ 
that can be compared with the WMAP results (Spergel et al. 2003; Peiris et 
al. 2003). 

Panel a) in Figure 5 presents the constraints in - dng/dlnk plane ob- 
tained from the analysis of the WMAP anisotropy measurements and the 
joint WMAP and Ly-a analysis from MCMC with dus/dlnk ^ 0. We found 
that both analyses favour a positive running of the slope dUg/dlnk ~ 0.023 
and an effective spectral index Ueff ~ —2.2 at A;p=0.05 Mpc~^. 
Panels b) and c) present the constraints in the Ues - crge"'^ plane and crge"'^ - 
dUes/dlnk plane obtained from MCMC with dUg/dlnk = ^. Prom the 
Monte Carlo Markov chain with dus/dlnk = we found Ues ~ —2.19 and 
dUes/dlnk 0.01. 

Our results differ from the value of the effective running of the slope dUg/ dink ~ 
-0.03 found by the WMAP team (Spergel et al. 2003; Peiris et al. 2003) at 

^ We have verified that the value of dUes/dlnk may depend quite critically on the 
choice of the step in k used to numerically compute the (central) derivatives of the 
power spectrum. For the considered models we find in practice quite stable results 
for steps in k less than ~ 5 % of fc, while using larger steps significantly affects the 
final results. 
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the same pivot wavenumber and indicate that a possible identification of a 
running of the slope, dng/dlnk ^ 0, at A;p=0.05Mpc~^ (multipole / ~ 700) 
with the current data is mainly an effect of the existing degeneracy in the 
amplitude-slope plane at this scale, the result being clearly consistent with 
the absence of running. 

In Figure 6 we compare the cosmological parameter constraints at 68% CL 
from the Ly-a analysis and the joint WMAP and Ly-a analysis with the 
cosmological parameter "simulated" constraints on cosmological parameters 
achievable by WMAP after 4 years of observations and by the combination of 
the three "cosmological" channels of Planck (Mandolesi et al. 1998, Puget et 
al. 1998, Tauber 2000) at 70, 100, and 143 GHz considering only the multipoles 
£ < 1500 and neglecting the Galactic and extragalactic foreground contami- 
nation. We assume Gaussian symmetric beams with the nominal resolution, 
a sky coverage of 80%, the cosmic variance and nominal noise sensitivity as 
sources of error, and neglect for simplicity possible systematic effects. Only the 
information from the temperature (TT) power spectrum is again considered. 
In the case of the "simulated" data we assume exactly the current WMAP 
data and error bars at i ^ 300, being the uncertainty in that multipole range 
dominated by the cosmic variance. We consider as fiducial model the best fit 
power law ACDM model to the WMAP data with dus/dlnk = (Table 1 
from Spergel at al. 2003). 

Note the improvement on power spectrum and reionization parameters achiev- 
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able by using the final WMAP data (improvement of about a factor of two) 
and that (of about a further factor of two) achievable with Planck by using 
only the temperature anisotropy data. Note also the role of Planck in re- 
ducing the error bar for the parameter w defining the equation of state of the 
dark energy component parameter. 

We find that adding the current Ly-cc information to the simulated WMAP 4- 
yr data only slightly reduces the error bars (of course, the relative improvement 
is significantly smaller by adding them to the simulated Planck data). 

4 Discussion and conclusions 

The recent detection of high values of the electron optical depth to the last 
scattering (Kogut et al. 2003, Spergel et al. 2003) by the WMAP sateUite 
(Bennett et al. 2003) implies the existence of an early epoch of reionization 
of the Universe at Zr ~ 20, fundamentally important for understanding the 
formation and evolution of the structures in the Universe. 
As the reionization is assumed to be caused by the ionizing photons produced 
during the early stages of star-forming galaxies and quasars, we evaluate this 
effect by computing the cumulative mass function of the high-z bound ob- 
jects for a class of fiat cosmological models with cold dark matter plus cos- 
mological constant or dark energy with constant equation of state, encom- 
passing different reionization scenarios. Our fiducial cosmological model has 
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Fig. 6. Cosmological parameter constraints at 68% CL from the Ly-a and 
WMAP+Ly-a analysis compared with the cosmological parameter constraints 
achievable by WMAP after 4 years of observations and by Planck. For all panels 
the meaning of the symbols is the same as in panel a). See also the text. 
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Qbh'^ = 0.024 ±0.001, fi^/i^ = 0.14 ±0.02, h = 0.72 ±0.05 as indicated by the 
best fit power law ACDM model of WMAP data (Spergel et al. 2003). 

Assuming that the virilization takes place at the collapse time and a constant 
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baryon/dark matter ratio in collapsed objects, we compute the fraction of the 
mass residing in gravitationally bounded systems as a function of the redshift 
of collapse at each linear scale and of the virial mass-temperature relation. We 
evaluate the formation rate of bound objects at ^ = 2.72 and their cumulative 
mass function was compared with the cumulative mass function obtained from 
the Ly-a transmission power spectrum (Croft et al. 2002). 
Our method allows to study reionization models by varying the amplitude, 
spectrum, and epoch of the reionization and the cosmological parameters. 
In the same time, as the high-z bound objects are rare fluctuations of the 
overdensity field, the tail of the cumulative mass function is sensitive to the 
rms mass fluctuations within the filtering scale cr{Rf, z). 
We find that the analysis of the cumulative mass function of the Ly-a systems 
indicates a reionization redshift in agreement with the value found on the basis 
of the WMAP anisotropy measurements, setting constraints on the amplitude 
of the power spectrum, (jg, similar to those derived from the X-ray cluster 
temperature function. 

Our joint analysis of Ly-a cumulative mass function and WMAP anisotropy 
measurements shows that a possible identification of a running of the slope, 
dus/dlnk 7^ 0, at A;p=0.05Mpc^^ (multipole / ~ 700) is mainly an effect of the 
existing degeneracy in the amplitude-slope plane at this scale, the result being 
clearly consistent with the absence of running, the other constraints based on 
WMAP remaining substantially unchanged. 
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We also shown that, for the set of cosmological models studied in this work, 
the error bars on the considered parameters can be reduced by about a factor 
of two by using the final WMAP data. The temperature anisotropy data from 
the forthcoming Planck satellite will further improve the sensitivity on these 
parameters, by another factor of two, and also the rehability of these results 
thanks to the better foreground subtraction achievable with the wider fre- 
quency coverage and the improved sensitivity, resolution and systematic effect 
control. This information jointed with the great improvement on the study 
of the Ly-a forest trasmission power spectrum (Seljak et al. 2002) achievable 
by the increase of the number of quasar spectrum measures expected from 
the Sloan Digital Sky Survey will allow to significantly better constrain the 
properties of the primordial density field at small scales. 
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